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Abstract 

We  compute  the  equilibrium  configuration  of  a  collision- 
free  plasma  contained  in  an  axially  symmetric  magnetic  field. 
The  plasma  is  characterized  by  a  non-scalar  pressure  tensor 
which  is  obtained  from  a  microscopic  distribution  function 
in  a  form  suggested  by  the  guiding  center  approximation.  The 
solution  is  calculated  in  the  limit  where  the  ratio  of  the 
width  to  the  length  of  the  plasma  region  and  the  ratio  of  the 
gas  to  the  magnetic  pressure  are  both  small.  Boundary  values 
at  the  midplane  as  well  as  the  shape  of  the  plasma  appear 
as  arbitrary  parameters  in  the  solution.  We  give  the  solu- 
tion to  a  corresponding  scalar  pressure  problem  for  comparison. 
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EQUILIBRIUM  CONFIGURATION  OF  A  PLASMA 
IN  THE  GUIDING  CENTER  LIMIT 

1.   Introduction 

We  present  a  calculation  of  the  equilibrium  configuration 
of  a  plasma  in  a  magnetic  field  for  a  plasma  with  a  non-scalar 
pressure  tensor.   Certain  other  studies  of  equilibrium  con- 
sider scalar  pressure  plasmas  (for  example  reference  [1])  as 
is  appropriate  for  a  collision  dominated  fluid.   Yet  whenever 
it  is  appropriate  to  characterize  the  plasma  as  a  collision- 
free  gas,  we  expect  that  the  presence  of  a  magnetic  field 
can  be  related  to  an  anisotropy  in  the  pressure  tensor.  For 
such  a  plasma  the  anisotropy  will  persist  since  there  are  no 
collisions  to  disturb  it. 

In  Section  2  we  assume  axial  symmetry  and  write  the 
equilibrium  equations  in  a  form  which  takes  advantage  of  this 
symmetry.  We  have  an  additional  symmetry  which  follows  from 
the  guiding  center  approximation:  the  pressure  tensor  is 
axial ly  symmetric  about  the  magnetic  field  lines.   That  is, 
in  a  coordinate  system  with  the  z-axis  tangent  to  the  magnetic 
line,  the  pressure  tensor  P  is  diagonal  and  has  the  form 
(see  reference  [2]) 


Px   0    0 


(1)       P  =   0    Pj_   0 

OOP 


■  / 
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where  Px  and  PM  are  called  the  components  of  P  perpendicular 
and  parallel  to  the  direction  of  the  magnetic  line.   In  general 
Px  is  not  equal  to  PM. 

The  equilibrium  equations  are: 


(2) 


V.  P  =  —  (V  x  B)x  B 
^o 


7  •  B  =  0 


v 


where  P  is  given  by  equation  (l),  B  is  the  magnetic  field, 
and  u_   is  the  vacuum  permeability.   We  need  an  additional 
relation  to  make  this  system  of  equations  mathematically  determ- 
inate since  there  are  four  equations  (2)  and  five  unknowns:  ?L, 
PM  and  the  three  components  of  S.   In  Section  3  the  additional 
relation,  similar  to  an  equation  of  state,  is  gotten  from  a 
microscopic  distribution  function.   The  general  form  of  this 
distribution  function  is  derived  from  the  guiding  center  approxi- 
mation.  We  ensure  the  validity  of  the  approximation  by  requiring 
that  the  magnetic  field  should  vary  slowly  with  position.   In 
addition  we  choose  a  specific  form  of  distribution  function 
which  is  useful  for  mirror  geometries.* 


* 

R.  F.  Post  (reference  [}])  has  reported  a  calculation  which 
somewhat  similar  to  ours.   He  does  not  use  a  distribution 
'unction  but  instead  chooses  a  form  for  the  magnetic  field  which 
approximates  some  experimental  cases.   He  then  solves  for  P,. 
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Section  4  contains  our  solution  to  the  equilibrium 
equations  (2).   We  make  two  simplifying  assumptions  which 
enable  us  to  exhibit  the  solution.  The  first  assumption  is 
that  the  ratio  of  the  width  to  the  length  of  the  plasma 
region  (the  aspect  ratio)  is  much  smaller  than  unity.   This 
assumption  serves  to  remove  the  magnetic  field  line  curva- 
ture from  consideration  and  to  reduce  a  partial  differential 
equation  to  an  ordinary  differential  equation.  The  second 
assumption  is  that  the  ratio  of  the  plasma  pressure  to  the 
magnetic  pressure  is  much  less  than  unity.  This  means  that 
the  presence  of  the  plasma  is  a  small  perturbation  on  the 
vacuum  solution  for  the  magnetic  field.  We  may  then  use  an 
iterative  technique  to  solve  the  equations.   We  find  that 
if  we  prescribe  a  linear  relation  between  Px  and  the  mag- 
netic flux  at  the  midplane,  as  well  as  several  other  standard 
boundary  conditions,  then  we  get  an  explicit  solution  to  the 
equilibrium  equations.   The  other  boundary  conditions 
(or  prescribed  values)  are:   the  prescription  of  the  plasma 
boundary  itself,  the  value  of  Px  at  the  origin  and  of  B 
at  the  boundary  of  the  midplane.  These  three  conditions 
(or  values),  subject  to  the  two  assumptions  above,  are  other- 
wise arbitrary  parameters  in  the  solution.  The  final  solution 
is  then  compared  with  the  solution  to  the  corresponding  scalar 
pressure  problem. 


-  6  - 


2.   Equilibrium  Equations  for  Axial  Symmetry. 

Since  we  deal  only  with  axially  symmetric  configurations, 
it  is  useful  to  introduce  the  cylindrical  coordinates   (r,  6,    z). 
(See  Figure  1  for  this  discussion. ) 

Fluxtubes,  V  =  constant 


axis 


Figure  1. 

We  assume  that  there  are  no  currents  in  the  z-direction 
in  the  equilibrium  state.   This  implies  that  Be  vanishes 
throughout  so  that  it  is  possible  to  introduce  a  stream  func- 
tion from  which,  by  differentiation,  we  obtain  B   and  B  . 

r  z 

Such  a   stream  function  is  the  magnetic   flux     i     defined  by 


(3) 


V(r,z)   =  2tt 


j 
o 


PdP  Bz(p,z) 


In  cylindrical   coordinates,   the   second  of  equations    (2)   is 
given  by 


CO 


V    •    B  = 


=  FcI^V+;|bz  =  0- 
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Prom  (2)  and  (4)  it  follows  that 
(5)  Br(r,z)  -  -  £;  H  • 


and 


(6)  V'-'l-K? 


From  (5)  and  (6)  we  obtain 


n    2 

(7)  (V  *  B)x  B  =  -(JU  (Ltf)frtf) 


where 


(8)  L^  =  — ^  -  F  ^  +  — ? 

dr  dz 


We  note  that   (VxB)x  B  lies  in  the  direction  of 
Slit.      It  seems  natural  then  to  evaluate  V  •  P  in  terms  of 
coordinates  along  <jf     and  B.   The  B-lines  lie  on  surfaces 
of  constant  f     (flux  tubes)  so  that  our  coordinate  system 
actually  consists  of  a  grid  of  the  surfaces  of  constant  tf 
and  the  orthogonal  trajectories  to  these  surfaces.  We  have 
already  assumed  that  P  has  the  form  given  by  equation  (1) 
In  any  other  rectangular  coordinate  system  (x^,  xg,  x^), 
the  components  of  P  are  given  by 


(9)  Pi(j  =  (P„  -  PJVj  +  ?L\y 
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where 


*-  B 


1*1   ' 

and 

fl,  if  i  =  J 

5i1  = 

J     0,  otherwise. 

We  shall  use  the  coordinate  s,  where  s  is  length  along 
a  B-llne,  and  write 


.   d  _  d 


and 


§=X(s)S 


where  >^(s)   is  the  curvature  of  the  B-line,  and  n  =W/|W 
is  the  principal  normal  to  the  B-line. 
We  get 


(10) 


V-  P=Us-(P»  "Px)-(P„  "Px)^|^}a 


+  (p„  "  Pi)  *U)n  +  VPj.  , 


where  B  =  |b| 
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Inserting  (7)  and  (10)  into  the  first  of  equations  (2) 
we  obtain  the  two  equations: 


(")        H*  -  (p.  -  O^P  -  o 


and 


(12,  *Ill  (P„  -  Px)  +  W  •  VP+ M^_ 

IV*  I  IV*|2      (2«)2n0 


We  can  write  equation  (11)  in  the  form 


(15>  d  log  B  =  P"  "  PJ-  ' 


where  the  indicated  differentiation  is  along  a  field  line. 

We  have  to  solve  the  two  equations,  (11)  and  (12),  for 
the  three  unknowns:   P„ ,  Px  and  if.      This  system  of  equations 
is  obviously  not  determinate  so  that  we  must  find  an 
additional  relation  connecting  the  three  unknowns.   We  do 
this  in  the  following  section  by  considering  the  individual 
particle  motion  in  the  guiding  center  limit. 
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3.   The  Form  of  the  Pressure  Tensor  In  the  Guiding  Center  Limit. 

In  this  section  we  find  the  additional  relation  linking 
Px,  P„  and  i     by  considering  the  individual  particle  motion 
of  the  ions  and  electrons  in  the  guiding  center  limit.   This 
limit  is  that  approximation  (3ee  reference  [4])  in  'which  the 
magnetic  moment  and  the  kinetic  energy  of  a  particle  are 
constants  of  the  motion  as  we  follow  the  guiding  center  of 
the  particle  along  a  magnetic  line.   This  means  that  if     is 
also  a  constant  of  the  motion  since  the  magnetic  lines  lie  on 
flux  tubes  (surfaces  of  constant  f).      The  guiding  center 
approximation  follows  from  the  restriction  that  the  Larmor 
radius  be  smaller  than  a  characteristic  length  of  the  system. 
We  assume  that  B  varies  slowly  with  position  so  that  the 
characteristic  length  is  the  maximum  of  (B/p?B|). 

The  following  arguments  apply  to  ions  and  electrons 
alike  so  that  subscripts  distinguishing  the  two  species  of 
particles  are  unnecessary.   The  equilibrium  equations  (2) 
combine  the  contributions  of  both  species. 

We  assume  the  existence  of  a  distribution  function  P 
for  the  particle  density.   The  distribution  may  be  written  as 
a  function  of  the  constants  of  the  motion.   In  the  guiding 
center  limit  this  motion  is  along  the  magnetic  lines  so  that 
it  is  convenient  to  erect  a  local  rectangular  coordinate 
system  on  the  magnetic  line  with  the  z-axis  tangent  to  B.   In 
this  notation,  the  particle  velocity  has  a  component  v„  along 
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the  z-axis  and  two  components,  designated  by  v\  ,  in  the 
plane  normal  to  B.  The  three  constants  of  the  motion  for 
a  guiding  center  particle  of  mass  m  are  the  energy  E, 
the  magnetic  moment  p,  and  the  flux  if     where 

(1*)  E=|(v2+V2), 

2 
(15)  n  =  ^jp  , 

and  if     identifies  the  magnetic  line  and  is  a  function  of 
position.  The  distribution  function,  F(E,  p,  it),     vanishes 
when  E  <  p£  since,  from  (14)  and  (15),  E  -  pB  =  5  v„  >  0. 
So  far  F(E,  p,  VO  applies  to  any  field  configuration. 
We  are  actually  interested  in  the  case  in  which  the  plasma 
is  confined  by  magnetic  mirrors.  The  distribution  will 
vanish  for  those  values  of  E  and  p  which  permit  the 
particles  to  escape  through  the  mirrors.  That  is:  let 
B   be  the  maximum  value  of  B  on  a  magnetic  line,  then 
particles  with  energies  E  >  pB   will  escape  from  the 
plasma  region.   In  summary  F(E,  p.,  if)   =  0  unless 


(E/Bm)  <  p  <  (E/B) 


Since  there  is  no  mass  motion  in  our  equilibrium  con- 
figuration we  can  write  the  components  of  the  pressure 
tensor  in  terms  of  F  as  follows: 
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(16) 


(17) 


P  {t)   =  m 


Px(f)  =| 


v^  F  dV, 


vf  F  dV, 


where  dV  indicates  integration  over  the  three  dimen- 
sional velocity  space.   In  terms  of  E  and  \i     these 
integrals  become: 

*    E/B  1/2 

(18)      P„(tf)  -  2cB   dE  J    dnF(E,  n,  tf)(E  -  |iB) 


E/B, 


m 


E/B 


(19) 


Px(^)  =  cB2   dE     nd[xF(E,  |1,  tf)(E  u  MB) 


-1/2 


E/B, 


rn 


where  c  =  2Tr(2m-5)~1' 2.   Following  a  particle  path,  we 
compute  from  (18)  and  (19)  that 


(20) 


Pn     —  p      p 

d  log  B  "  r"  '   •>- 


Equation  (20)  is  the  same  as  equation  (15)  although  the 
latter  was  derived  from  a  macroscopic  momentum  balance. 
We  remark  that  equation  (20)  does  not  depend  on  our 
restriction  to  mirror  geometries.   It  is  a  general  con- 
sequence of  the  guiding  center  approximation.   Had  we 
chosen  to  derive  equation  (20)  before  requiring  that  F 

vanish  for  |1B  <  E,   the  result  would  have  been  the 

m 
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same  since  the  lower  limit  in  the  i-i-integration  in  (18) 

and  (19)  would  have  been  zero  instead  of  E/B  . 

m 

In  general,  the  specification  of  F  may  be  quite 
arbitrary.   In  any  particular  instance  it  may  depend  in  a 
complicated  way  upon  design  parameters  such  as  the  manner 
of  injection.   One  could,  in  principle,  inject  an  arbitrary 
number  of  particles  with  given  u-  and  E  on  each  magnetic 
line.   Each  arbitrary  specification  of  P  determines  P„ 
and  Pj_  as  functions  of  ii .      Among  possible  representa- 
tions of  P  we  might  naturally  consider  terms  of  the  form 
f(E,  ^)(u-B  -  E)n,  which  fall  to  zero  properly  as  E 
approaches  M-B  ,  or  even  a  sum  of  such  terms: 


00 
P(E,  M-,  1>)    =  Y~  fn(E,*)(HBm  -  E)n. 


In  order  to  obtain  a  fully  worked  out  example  we  shall 
specify 

^(E^MHB,*-  E)n, 


m 


(21)      F(E,  M-,  if)    =  < 


for  (E/Bm)<  H  <  (E/B), 


0,  otherwise, 


where  n  >  0  is  not  necessarily  an  integer. 

With  this  choice,  equations  (18)  and  (19)  become: 

Bm  n        (n+3/2) 

(22)  P„(*)  =  A(*)(^)  (1  "  I") 

m 
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B  n     „   (n+l/2)  R 

(25)  Px(*)  =  A(*)(g2i)  (1  -  I;)       ((n+D  +  ^ 


m  m 


where 


A(*)  _  2c  r  (3/2)  R"+l)  f  f(E,*)E(n+5/2)  dE 
Hn  +  5/2)     J 


* 


and   P(x)  denotes  the  gamma  function. 

In  (22)  and  (23),  only  A(tf)   depends  on  whether  we  have 
electrons  or  ions.   We  add  the  contributions  of  the  two  sets 
of  particles  to  form  a  new  A(tf)  and  assume  henceforth  that 
equations  (22)  and  (23)  refer  to  the  total  plasma. 

Let 

(24)  M(tf,  s)  =  Bm(l/0  /  B(tf,  s), 


then 


(25)  Pi(^,  s)  m   (n-»-l)M(*,  s)+  1/2 

P„(tf,  s)       M(tf,  s)  -  1 

(The  coordinate  s  will  be  omitted  henceforth.)  Equation  (25) 
is  the  relation  which  was  needed  to  make  the  system  of  equa- 
tions (2)  determinate.   From  (25)  we  see  that  (Pj/P,, )  >  (n+l). 


See  for  example  formulas  48l  and  482  on  Page  62  in  B.  0. 
Peirce,  A  Short  Table  of  Integrals,  Third  Revised  Edition, 
Ginn  and  Company. 
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In  solving  the  equilibrium  equations  we  will  prescribe 
?L     as  a  function  of  1/     at  the  midplane  s  =  0.  It  is 
convenient  to  give  A(V)  in  terms  of  Px     (the  subscript  o 
indicates  values  at  the  midplane)  and  obtain  equations  (22) 
and  (25)  in  a  more  useful  form.  From  (23)  we  have 

(26)      A(*)  =  Plo(^)M;n(l-l/Mo)-(n+1/2)[(n+l)+  1/(2M0)]-1, 
Substituting  for  A(^)  in  (22)  and  (23)  we  get 


(27)  Pl(*)   -  Plo(*)[«  )n6=$l 


(n+1/2) 


n+l)+  1/(2M) 


mm 


>+!)+  V(2Mp 


and 


(28) 


P,,(*)  = 


Pl0(^)(^)n(i-l/M)(n+5/2)(l-l/M0)-(n+1/2)[(n+l)+l/(2M0)]-1. 

From  these  equations  we  may  draw  a  few  qualitative  results. 
We  see  that  P„  and  PA  decrease  monotonically  as  B 
increases  and  finally  vanish  at  the  mirrors  where  B  is  a  maxi- 
mum.  PM  vanishes  faster  than  PA  at  the  mirrors.   If  Pj.0(VO 
should  vanish  for  some  V  =  T,     then  PA  =  PM  =  0  on  the 
entire  flux  tube  f   .  This  is  used  as  the  natural  way  of 
defining  the  boundary  of  the  plasma. 
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4.   Equilibrium  Configuration  for  Low  Aspect  Ratio  and  Low  p 

The  problem  given  by  the  equilibrium  equations,  (12)  and 
(13),  and  the  relation  (25)  is  mathematically  determinate 
once  suitable  boundary  conditions  are  given.   The  discussion 
in  the  last  section  showed  that  Plt  and  Pj_  already  satisfy 
equation  (13)-   This  was  a  consequence  of  the  guiding  center 
approximation.   We  will  solve  the  remaining  equation  (12) 
after  making  three  simplifying  assumptions.   The  first,  which 
is  not  very  restrictive,  is  that  the  system  is  symmetric  about 
the  midplane.   The  second  is  that  the  aspect  ratio  (the  ratio 
of  the  width  to  the  length  of  the  plasma  region)  is  much 
smaller  than  unity.   This  considerably  simplifies  equation 
(12).   The  third  assumption  is  that  (3   (a  representative  ratio 
of  gas  to  magnetic  pressure)  is  much  smaller  than  unity.   This 
means  that  we  can  treat  the  plasma  as  a  small  perturbation 
imposed  on  a  vacuum  magnetic  field.   Our  solution  which  will 
be  explicit  will  be  given  to  first  order  in  0. 

We  need  suitable  boundary  conditions  on  ii     before  we 
can  solve  the  equilibrium  equation.   First  we  introduce  the 
notation  to  be  used  in  the  rest  of  this  section.   Quantities 
at  the  midplane  (see  Figure  2)  are  given  the  subscript  o, 


£ 


->-  z 


zn> 


o  z 

\  m 

midplane 


Figure  2 
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those  at  the  mirrors  the  subscript  m.  Quantities  at  the 
plasma  boundary  will  have  a  bar  over  them,  such  as  J     or 
B.   In  order  to  develop  the  boundary  conditions  we  first 
prescribe  the  plasma  boundary  r(z).   The  only  requirements 
on  r(z)  are  that  it  vary  slowly  with  z,  that  it  have 
its  maximum  at  the  midplane   (z=0),  and  that  it  decrease 

to  a  minimum  at  the  mirrors   (z  =  *  z  ).   The  plasma 

m       ^ 

boundary  can  be  related  to  the  coil  system  which  is  at  our 
disposal.  We  next  prescribe  a  linear  relation  between  ?x 
and  y     at  the  midplane : 

(29)  pi0(V  = a  (vF<0))-  Vr)]  ' 

where  a  and  y   (r(o))  are  constants  at  our  disposal. 
In  addition  we  give  B   (B  at  the  boundary  on  the  mid- 

plane)  and  PXo(0)   (pio  at  the  oriSin). 

The  choice  (29)  leads  to  an  explicit  solution  at  the 
midplane.   Then  we  determine  J   =  y   (r(0)).   This  is  one 
of  the  boundary  conditions  on  y(r,z).      The  other  comes 
from  equation  (3)  which  implies  that  y     vanishes  on  the 
axis.   In  summary,  the  boundary  conditions  on  y(r,z)     are: 


(30) 


y{o,z)  =  o 

y(r{z),z)   =  ? 


<. 
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From  (29)  and  (30)  we  see  that  Pio(^)  vanishes  on  r(z). 
This  in  turn  implies  that  Px  and  Pn,   given  by  (27)  and 
(28),  also  vanish  on  the  plasma  boundary.   Instead  of  pre- 
scribing  r(z)  we  could  just  as  well  have  prescribed  B  on 
the  axis  as  an  increasing  function  of  z  with  maxima  at  the 
mirrors.   We  would  then  give  the  quantities  r  (0 )   and 
Px  (0)  and  employ  the  new  boundary  conditions: 


r 


(30') 


TH0,z)  =  o 


JL|i  |     =  B(0,z) 
27TT  dr    '   q     zv  '  ' 


The  second  equation  in  (30' )  follows  from  (6)  and  the  fact 

that  B  vanishes  on  the  axis.  The  two  different  sets  of 
r 

boundary  conditions  yield  substantially  the  same  solutions. 

From  (24)  and  (27)  it  is  evident  that  we  must  first  find 
B  (1>) ,   the  maximum  value  of  B  attained  at  the  mirrors, 
before  we  can  solve  in  the  interior.   It  turns  out  that  B 
is  a  constant  for  all  of  the  magnetic  lines  since  the  plasma 
vanishes  at  the  mirrors. 

We  solve  the  equilibrium  equation  at  arbitrary  z  to 
first  order  in  3  by  a  simple  iteration  process  to  be  des- 
cribed at  the  appropriate  place.   This  solution  is  then  compared 
with  the  solution  to  the  corresponding  scalar  pressure  problem. 

The  detailed  calculations  now  follow. 
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Low  aspect  ratio 

We  define  the  aspect  ratio  €  by 


(31)  €  =  F(o)/zm 


and  assume  that  e  «  1.  This  implies  that  all  functions  of 
interest  may  be  written  in  the  form 

f  =  f(r,  ez). 

p 
We  neglect  all  terms  of  order  €   such  as  the  field  line 

curvature   ,   second  derivatives  with  respect  to  z,  and 

products  of  first  derivatives  with  respect  to  z.   To  this 

approximation  equation  (12)  becomes 


/,«x  6^    1  dtf     u    2    2  dPx  /  dtf 

(32)  ^2"  "  F  dT  "  "47r  ^or  Sr-^SF 


We  note  that  z  enters  only  as  a  parameter  in  the  differen- 
tial equation.   The  formal  solution  to  an  equation  of  the 
form 


(33)         4-fS=^' 
dr^   r  dr  . 


is 

r     t 
(34)  tf  =  4tt2   tdt  I  sg(s)ds  +  Cr2, 

0     0 
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where  an  additive  constant  has  been  set  equal  to  zero  since 
f     vanishes  on  the  axis.   From  our  assumption  that  £  is 

also  much  smaller  than  unity  we  shall  show  that  g(r)  will 
be  independent  of  r.   For  this  case  (3^)  reduces  to 


4 
(24-)         1/   =  7r2g  |-  +  Cr2. 


The  constants  will  be  determined  from  the  boundary  conditions 
as  functions  of  z. 

A  further  consequence  of  the  assumption  that   e  be  small 

is  that  B  is  equal  to  B   to  order  €.   This  follows  from 

o 
(5)  since  B   is  of  order  e     but  appears  only  to  order  e 

in  the  magnitude  of  B.  We  evaluate  B  by  equation  (35) , 

henceforth: 


C5)  B(r,,).^g 


Solution  at  the  midplane 

To  find  the  solution  at  the  midplane  we  substitute  (29) 
into  the  right  hand  side  of  (32)  and  obtain 

dV    -,    dtf      ,, 
(36)  ^-J^-^V 

From  equation  (3^')  the  solution  is 
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4 
(37)  *Q(r)  =  ir\a   f-  +  Cr2. 


We  prescribe  PlQ(0)   (p1q  at  the  origin)  and  BQ  (B  on  the 
boundary  of  the  midplane).  The  constants  C  and  a  are 
determined  from  (29)  and  (35). 


(38)  C  =  ttBo(0), 


and 

B  (  0)        .  . 
(39)  a  =-44-  t(l+P)l/2-  U* 

^oro 


where 

(^0)  0  =  2^oPio(°)  /  Bo(o)> 

(^D  BQ(0)  =  (B2  -  2noPio(0))l/2, 


and 


rQ  =  r(0) 


Low  g  limit 

The  quantity  p,   defined  by  (40),  is  the  maximum  ratio 

of  gas  to  magnetic  pressure  for  the  configuration  we  have 

2 
chosen.  We  assume  that  6  «  1,  ignore  terms  of  order  p  , 

and  obtain: 

-  22  - 


fB.(O) 

(59')         a  =p_^|p, 


(42)  *Q(r)  =  7rBo(0)r2(l  +||«-)  | 

and 

(43)  ?  =  ^0(Fo)  =  ttBo(0)Fo2|i  +  £j. 

The  boundary  value  ?  given  by  (43)  is  used  in  (JO)  to  deter- 
mine the  solutions  off  the  midplane. 

For  future  use,  we  obtain  from  (35) 


(»)  B0(r)  -Bo(0)jl  +  (jrf  f 


Solution  at  the  mirrors. 

We  need  B   (=  B(r,z_))  in  order  to  evaluate  P.   for 
m   x    v   m  *■ 

z  ^   0.   At  the  mirrors  P  vanishes  so  that  we  must  solve 

(32)  with  the  right  hand  side  set  equal  to  zero  and  impose 

the  boundary  conditions  (30).   We  then  determine  B  from 
(35)  and  obtain 


(^5)  Bm  =  7/(0 


m    '  v   m 


where 


r  =  r(zm) 

m    v  m' 
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The  quantity  Bm  is  a  constant  for  all  flux  tubes  to 
this  approximation. 

Solution  off  the  midplane. 

The  right  hand  side  of  equation  (32)  is  only  of 
order  B.   Therefore  it  is  reasonable  to  use  an  itera- 
tion scheme  to  solve  this  equation.  The  first  step  is 
to  obtain  the  solution  with  the  right  hand  side  of  (32) 
equal  to  zero  and  with  the  boundary  conditions  (30). 
From  this  solution,  denoted  by  i^   (r,z),  we  find 
B^°'(r,z)  by  means  of  (35)-   The  next  step  is  to  insert 
tT°)  and  B^0)  into  the  right  hand  side  of  (32)  and 


solve  for  the  next  approximation  to  V  which  is  denoted 

cti< 
,(D 


by  ^  '(r,z).   Further  corrections  will  be  of  higher  order 


in  (3  so  that  we  stop  at  ijf 

We  need  equation  (27)  in  the  more  convenient  form 

(46)  Pj.(r,z)  =  PiQ(tf)  •  h(r,z), 

where 


(47)      h(r,z)  = 


'bq(*)  ^ 

B(r,z)y 


nr  \(n+l/2), 

/l-[B(r,z)/Bm]        j (n+l)+[B(r,z)/2Bm] 


^l-[B0(^)/Bm]J      l(n+l)+[B0(^)/2Bm] 
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The  notation  B  (VO  means  that  B   is  evaluated  at  the  mid- 
o v  '  o 

plane  on  the  same  flux  tube  1/     as  the  point   (r,z).   The 

function  h(r,z)   is  unity  at  the  midplane  and  decreases  with 

increasing  B  to  zero  at  the  mirrors.   We  will  find  that  the 

terms  in  r  in  B(r,z)   (just  as  for  B  (r))   are  only  of 

order  B«   The  term  on  the  right  hand  side  of  (32)  is  already 

of  order  B  so  that  we  need  only  evaluate  h  on  the  axis. 

This  formal  mathematical  procedure  is  reasonable  whenever 

the  mirror  ratio  B  /B   is  not  very  close  to  unity.   We 

mo 

insert  B^0^  and  t^°)   into  (47)  and  get  h^°^(0,z).   The 
equation  for  ^  '(r,z)   is 

(48)      iV^-I^=4.Vuoah(°>(0,z). 

In  order  to  find  i/^° '  and  V^  we  need  only  the  formal 
solution  (3^')  and  the  boundary  conditions  (50).  Equation 
(35)  then  gives  B*°^  and  B^  '.   The  results  are: 


(49)      *(0)(r,z)  =  j(^jjf    , 


(50)      B(o)(r,z)  =  B(o)(z)  =  J/(m   (z)) 


1 


/tr/-\i2^n+sV„,  .„i.r=  f=!     m2 


(5,  h(°)(o,zhfeM  W!l'  a'"1'^'','' 

ro  /        i-r? vf  i2  ?rn+iurF  /f  i2 


^lAc/J  2(n+l)  +  [rm/?of 


'  > 
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(52)     tf(1)(r,z)=^°>(r,z) 


2   > 


i  +  £h<°><o,2)  m       '   .  : 


r0   |\r(z) 


(53) 


B(1)(r,z)=B(o>(Z/ 


l+|h(°)(0,z)(!l£lj|2 


/ 


m] 


We  could  use  (42),  (44),  (52)  and  (52)  to  evaluate 
PL     numerically  as  a  function  of  r  and  z.  To  lowest 
order  in  p,  however,  it  is  sufficient  to  use  1^'°' 
and  h(o)(o,z)  in  (46).  Using  (29),  (39'),  (43),  (46), 
(49)  and  (51)  we  find  that 


"\ 


(54)       Mr  z)   =  Ji  -I 


(B^(0)/2n0) 


ir(z) 


h(o)(0,z), 


and  from  (25)  that 


ri-(rm/r(z))2 


1-1 


(55)     -^^^L—  =  2^<   

(B2(0)/2no)     [2(n+l)+(rm/r(z))^ 


•h 


(o)(0,z) 


Equations  (52),  (53),  (54)  and  (55)  (where  h^  is 
given  in  (51))  constitute  our  solution  to  the  equilibrium 
equations  to  order  p.   The  plasma  boundary  r(z)  as 
well  as  P.  (0)  and  B   are  arbitrary  parameters  in  the 
solution,  subject  to  the  restrictions  that  e  «  1  and 
that  0  «  1. 
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The  corresponding  scalar  pressure  problem. 

In  order  that  the  scalar  pressure  problem  correspond 
to  the  previous  problem  we  preserve  the  same  values  at  the 
midplane,  the  same  plasma  boundary  r(z),  and  the  same 
linear  relation  (29)  between  the  scalar  pressure  Pg  and 
tf  at  the  midplane.   From  equation  (15)  it  follows  that  Pg 
is  constant  on  a  flux  tube  since  P„  -  Px.   We  have  therefore 


(56) 


Ps  =  a{J-i/{r,z)) 


throughout  the  plasma  region.   Equation  (32)  becomes 


(57) 


d2*, 
i? 


M 


-h&-  ^£v 


where  the  subscript  s  denotes  the  scalar  pressure  quanti- 
ties.  We  haye  the  same  boundary  conditions  (30)  for  this 
problem.   The  solution  to  (57)  is  given  by  (3V).   We  obtain 


(58) 


(59) 


*«  =  ? 


\ 


B_  = 


r(z) 


J 


i+e  £M  L 


>(z 


V 


- 1 


7*2(Z) 


l*«(¥ 


2!ifo 


-  1 


These  solutions  hold  for  arbitrary  £  although  for  the  com- 
parison we  again  take  £  «  1.   To  the  same  order  of 
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approximation  as  before,   the   scalar  pressure  is  given  by 


(60)  — J5— § =  p 

(BJj(0)/2n0) 


/■ 


1   - 


r(z) 


The  comparison  between   (V^  ' ,   B^  ')  and   (^  ,  B  )   is 

s   s 

not  particularly  revealing  since  they  differ  by  a  term 
of  order  B.   In  fact  by  setting  h^o'(0,z)  equal  to 
unity  in  (52)  and  (53)  we  get  the  scalar  pressure  results 
(58)  and  (59).   The  function  h^°'  contains  the  essential 
difference  between  the  scalar  pressure  and  non-scalar 
pressure  results. 

Perhaps  the  more  interesting  comparison  is  the  ratio 
R  =  Pj/P  .   From  (54)  and  (60)  we  obtain 


(61)  R  =  h(o)(0,z) 


This  ratio  starts  out  equal  to  unity  at  the  midplane  and 

decreases  monotonically  with  increasing  B  until  we  reach 

the  mirrors  where  R  vanishes.  This  is  to  be  expected 

since  P  remains  constant  on  the  flux  tubes  while  P 
s 

and  PL  both  vanish  at  the  mirrors. 
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